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Abstract 

Recent work on the complete wetting transition for three dimensional systems with short- 
ranged forces has emphasized the role played by the coupling of order-parameter fluctuations 
near the wall and depinning interface. It has been proposed that an effective two-field 
Hamiltonian, which predicts a renormalisation of the wetting parameter, could explain the 
controversy between RG analysis of the capillary-wave model and Monte Carlo simulations 
on the Ising model. In this letter results of extensive Monte Carlo simulations of the two-field 
model are presented. The results are in agreement with prediction of a renormalized wetting 
parameter ui. 



There are a number of long standing controversies in the study of phase equilibria at fluid 
interfaces, related to wetting transitions for systems with short-ranged forces at the marginal 
dimension d = 3 (for a recent review see for example Q). The standard model used to describe 
fluctuation effects at wetting transitions is the effective interfacial Hamiltonian (also called the 
capillary wave (CW) model) 

•£ Q/3 (T) 



where ^(y) > is a collective coordinate which represents the distance of the fluctuating fluid 
interface (separating bulk phases a and (3, say) from a wall situated in the I = plane. Here 
X Q/ 3 is the interfacial stiffness coefficient and W(£) is the effective binding potential which is 
usually specified as 

W(£) = Ae~ Kl + Be- 2K£ + hi (2) 

where k = 1/^ is the inverse bulk correlation length and h is proportional to a bulk ordering 
field. At a critical or complete wetting transition the mean interface displacement (I) (and 
other lengthscales) diverges as the external field approach critical values. For critical wetting, 
B is positive and the transition occurs at h = A = 0, corresponding to T — > T\y, the wetting 
temperature, while for complete wetting A > 0, i.e. T > Tyy and the transition occurs for h — > 
for Tc > T > Tw- 

One well known controversy concerns the critical wetting transition. For this case, RG calcula- 
tions based on (|l]) (see for example [Q] ) , predict strong non- universality for the critical exponents 
and amplitudes depending on the value of the wetting parameter ui{T) = ksTn 2 j ' \.TxYi a p at the 
transition temperature Tyy. However, extensive Monte Carlo simulations of wetting in a three 
dimensional Ising model with nearest-neighbour interactions by Binder, Landau and Kroll ||, 
show only small deviations from mean-field theory consistent with uj ~ 0. This contrasts sharply 
with the expected value u ~ 0.8 (at the temperature where the simulations were performed). 
Prompted by the discrepancy Gompper and Kroll [GK] |5j performed very accurate Monte Carlo 
simulations on a lattice version of (|l]) and found very good agreement with the renormalisation 
group treatment, leaving the controversy open. 

More recently attention has focused on the complete wetting transition. For this case CW theory 
predicts that, while the critical exponents retain their mean-field values in the limit of h — > 0, 
critical amplitudes are renormalised. For example, the mean distance {£) diverges at leading 
order as 

k(£) « e\n(l/h) (3) 
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where 6 depends on oj. In particular for oj < 2, CW theory predicts 6 = 1 + ^. A new 
controversy appeared with the Ising model simulations of Binder, Landau and Ferrenberg 
(at oj ~ 0.8) consistent with 6 ~ 1.6 — 1.8, while the CW model predicts only ~ 1.4. 

These results serve to emphasize that wetting at the marginal dimension is a rather special 
case and that the coarse-grained Hamiltonian ([!]) may not capture all the essential physics. 
Note that no such problems arise in two dimensions where the CW model yields predictions for 
universal critical exponents, amplitudes and scaling functions which are in excellent agreement 
with the exact Ising model results ||. Motivated by this, several authors have undertaken a 
careful examination of the foundations of the RG theory and of the validity of the interfacial 
Hamiltonian. In particular Fisher and Jin [FJ] j?]] initiated a method of systematically deriving 
interfacial Hamiltonian from an underlying Landau-Ginzburg- Wilson [LGW] model. For critical 
wetting, they suggested that (||) should be modified by allowing a position dependent stiffness 
coefficient and as a result argue that the transition may in fact be very weakly first-order. More 
recently, Parry and Boulter [PB] || argued that this is still not a sufficient modification of the 
CW model, since it does not account for fluctuations of the order parameter near the wall and 
depinning interface. Following the methods by [FJ] , they introduced a two-field Hamiltonian 
for complete wetting 

H[h,h} = J dy{^ v (h,e 2 )Ve^V£ 1/ + U(h) + W(£ 2 -£ 1 )} (4) 

which is a functional of two collective coordinates ^i(y) and £ 2 (y) (representing the locations of 
generalised surfaces at the wall and the a(5 interface respectively) and 5] is a position dependent 
stiffness matrix ||. The potential U(£\) binds the lower surface to the wall and since the 
fluctuations of £\ are small can be approximated by U{£\) = |uo^i- Here v$ = Eu/^L, where 
£ w g is the the finite correlation length associated with the intrinsic fluctuations of the order 
parameter near the wall. The fluctuations of ^i(y) are not included in the CW and FJ models. 
The term W{£ 2 —£\) is similar to the binding potential in @. In the limit of complete wetting, 
£ 2 unbinds from the wall while £\ remains bound. In their RG analysis, [PB] show show that 
the coupling to a weakly fluctuating field results in an effective value of the wetting parameter. 
In the limit vq — > oo or £ w g — > the fluctuations of the lower surface are completely suppressed 
and the CW/FJ result 6 CW = 1 + % is recovered. In the other limit, vq — > or £ w p — > oo, the 
Hamiltonian (^]) can be written in terms of center of mass coordinates which gives an effective 
stiffness coefficient S" 1 = S^ 1 + and results in a capillary parameter oj e ff = oj 2 + oj\. For 



intermediate values of vq, linear Q and non- linear RG analysis [10] predict that the effective 
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capillary parameter is 



Consequently, the critical amplitude 9 is given by 9 = 9 CW + A9 with 



2V4jt(E 11 +o /A 2 ). 

where A is a momentum cutoff. They estimated || A9 ~ 0.3 ± 0.1, which gives the prediction 
9 ~ 1.7 ±0.1, in reasonable agreement with the Ising simulations in the complete wetting regime. 
For the subsequent development in theory, it seems important that [PB] theory is tested by 
independent methods. The crucial prediction is that the capillary parameter is renormalised 
due to the coupling to a field ^i(y) which only weakly fluctuates. To this end, in this letter, 
the Monte Carlo studies of critical wetting by [GK] || have been generalised to the two-field 
Hamiltonian in the complete wetting regime. 

The simulations were performed on a lattice version of (|j). The two surfaces l\ and £2 are 
modelled by two parallel Lx L square lattices with periodic boundary conditions. The variables 
x(i,j) and y(i,j) are the distances of the site of the lower (£±) and upper (£2) surface 

respectively, from the wall (located at zero). The uncoupled gradient terms are approximated 
by 

77 — J2 ~ x s(hj)) 2 , -77, — ^2 (y(i,j) -ys(i,j)) 2 (7) 

Wttuji f- / 167ru;2 , 

where 5 indicate the four nearest-neighbours of the site and the sum is over all sites 
for i,j = 1, L, and x and y are measured in units of the bulk correlation length £&. To test 
the program, the potential (||) is first chosen as A = 0,B = 1 (critical wetting), En is set very 
large to suppress the fluctuations of the lower surface i\ , the cross-gradient coupling term S12 is 
kept to zero and u>2 takes the values given by [GK] ||. The results are in very good agreement 
with their founding. For complete wetting, the potential is chosen as A = 1 and B = and 
the cross-gradient coupling term £12 kept to zero ||. The lattice size L is varied between 10 
and 40 to study finite-size effects. The SOS approximation is adopted and x(i,j) and y(i,j) 
are treated as continuous variables to avoid problems related to roughening. For simplicity we 
ignore the exponentially small position-dependence of the stiffness matrix elements since they 
play no role in determining the increment A.9 pj. The hard wall condition (y > x) is used here, 
although the prediction (|5|) is also valid in a soft wall approximation [||, (To]] (see for example j|] 
for a discussion regarding the CW model). The wetting parameters are chosen as uj\ = 1 and 
LO2 = 0.8, and so the maximum increment to 9 is A9 max = 0.5 (corresponding to vq = 0). 
A Metropolis algorithm has been used on a Dec Alpha 3600. At each Monte Carlo sweep, one 



of the two surfaces is chosen at random 2L 2 times, and then a randomly picked site is 
updated by a random number in the interval [— z m , +z m ]; z m is set so that approximately 50% 
of the updates are accepted (but other values have also been used to test the convergence rate). 
The first 10 4 sweeps are discarded to thermalise the system and 10 6 sweeps are used to calculate 
the averages. The block average method is used to obtain statistical errors and the integrated 
autocorrelation time r% n t is calculated as in Madras and Sokal |DJ . For example for the parameter 
vq = 100 and field h = 0.1, Ti n t = 15 and for h = 0.01, T{ n t = 113, while for vq = 0.05 and 
field h = 0.1, Ti n t = 25 and for h = 0.01, Ti n t = 159. For a given value of the parameter vq the 
average distance (y) of the upper surface from the wall, is calculated for several values of the 
fields h. Firstly, to approximate the limiting value vq « oo or « 0, vq is set at vq = 100. In 
Figure 1, (y) is plotted against ln(l//i) as in eq.(^). The values of (y) lie on a straight-line and 
a logarithmic regression gives a slope 9 = 1.43 db 0.008. The other limit, Vq ~ or ~ oo, 
is approximated by vq = 0.05. For smaller values of Vq the fluctuations become very large, the 
autocorrelations increase and the simulations are not reliable. The values of (y), plotted against 
ln(l//i) in Figure 1, show that they lie on a straight-line, but now with a slope 9 = 1.65 ± 0.02. 
Figure 1 suggest that 9 is vq dependent and other intermediate values of vq are chosen to show 
this. In Figure 2 are plotted the fitted values of the slopes 9 corresponding to different values 
of vq for lattice sizes L = 20, 40. Importantly, the graph shows that the value of 9 is indeed 
larger than the CW prediction 9 CW ~ 1.4 for finite values of vq, consistent with (|6]). Due to the 
dependence on the cutoff A a precise quantitative comparison between simulation and theory is 
not possible. Nevertheless a fit of the numerical results to the theoretical prediction @ yields 
a value for the cutoff A close to unity, which is encouraging. We also note that as t>o increases, 
the increment A# decreases and we recover the CW model prediction 9 CW ~ 1.4 as Vq — > oo. 
For small vq the measured value of 9 is considerably bigger than the CW result and increases 
non-linearly with 1/vq. As mentioned above, the limiting value should be 9 = 1.9, but this can 
not be reached in the simulations since it corresponds to infinite fluctuations. Unfortunately, 
extrapolation to v$ = is not possible either, given the non-linear dependence of 9 on vq and 
large finite-size effects which would occur in this limit. Nevertheless the observed increase of 9 
with vq testifies to the strong influence of the coupling between the two fields and is certainly 
consistent with the basic prediction of the coupled Hamiltonian theory. 

In the light of [PB] theory ||, the results are interpreted as follow: for very large values of 
vq the lower surface is very stiff and it does not fluctuate, hence the capillary wave result of 
qCW ~ i_4 j s recovered. For decreasing values of vq, even small fluctuations in the lower surface 
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coupled to the upper surface produce a change in the critical amplitude 6 as predicted by the 
theory and also in agreement with the Ising model simulations. This is the first time that 
RG predictions of an effective Hamiltonian, simulations on the Ising model and simulations of 
the same effective Hamiltonian are in qualitatively and quantitatively agreement. It would be 
interesting to extend this work to critical wetting and investigate if (Q) can shed light on the 
controversy there. However this is beyond the scope of the present letter. 
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Figure 1: Mean interface displacement (y) against log(l//i), for vo = 100 (circles) and vo = 0.05 
(squares). The error bars are within the symbols. The solid lines are the logarithmic regressions. 
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Figure 2: Critical amplitude 9 as a function of the inverse parameter vo, for lattice size L = 20 
(circles) and L = 40 (squares). 
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